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Numerical study of the one-dimensional Frenkel Hamiltonian with on-site randomness is carried 
out. We focus on the statistics of the energy levels near the lower exciton band edge, i. e. those 
determining optical response. We found that the distribution of the energy spacing between the 
states that are well localized at the same segment is characterized by non-zero mean, i.e. these states 
undergo repulsion. This repulsion results in a local discrete energy structure of a localized Frenkel 
exciton. On the contrary, the energy spacing distribution for weakly overlapping local ground states 
(the states with no nodes within their localization segments) that are localized at different segments 
has zero mean and shows almost no repulsion. The typical width of the latter distribution is of the 
same order as the typical spacing in the local discrete energy structure, so that this local structure is 
hidden; it does not reveal itself neither in the density of states nor in the linear absorption spectra. 
However, this structure affects the two-exciton transitions involving the states of the same segment 
and can be observed by the pump-probe spectroscopy. We analyze also the disorder degree scaling 
of the first and second momenta of the distributions. 



I. INTRODUCTION 

Since the pioneering work by Mott and TwoseEl it is 
well known that all states in one dimension (ID) 
localized in the presence of uncorrelated disorder 
which means that a quasi-particle wave function has a fi- 
nite amplitude within a finite space interval and vanishes 
outside. The size of this interval (localization length) in- 
creases with the decrease of disorder. The ID localiza- 
tion theorem was supported latei^by the one-parameter 
scaling theory of Abrahams et aln (see a comprehensive 
review in Ref. 

The concept of ID localization have been successfully 
applied to explanation of the optical properties of linear 
molecular aggregates and conjugated polymers, in which 
the elementary excitations are Frenkel excitons (for a 
comprehensive review see Refs. |,|,|). The one-exciton 
states that are closer to the bottom of the exciton energy 
band couple better to the light and hence determine the 
linear optical response. For this reason the energy (and 
wave function) statistics of the lower states has been at- 
tracting great deal of attention. In Ref. |l^ the heuristic 
concept of the hidden low-energy structure of a ID local- 
ized exciton was put forward. This concent .aias-.partly 
supported by the numerical simulations .EJIIj'Ejo How- 
ever, it has not been proved directly by detailed statisti- 
cal study. 

The basic idea of the hidden low-energy structure con- 
cept is as follows. The lower energy one-exciton eigen- 
functions obtained for a fixed realization of disorder can 
be grouped in sets of two (or sometimes more) states 
which are localized at the same segment of the linear 
chain. These segments are almost non-overlapping. The 
typical length of a segment N* depends on the degree of 



disorder (the quantity TV* is, often called the number of 
coherently bound moleculeala). The energies of the states 
belonging to the same set are well separated, so that they 
form a local discrete energy structure. The mean energy 
of such local manifolds is distributed within the energy 
interval of the same order as the typical energy spacing 
between the levels of a local manifold. For this reason 
the local energy structure appears to be hidden; it does 
not manifest itself neither in the density of states nor 
in the linear absorption spectra. Nevertheless, it can be 
revealed frqnj, tbe. optical response involving two-exciton 



transitions .Eilii 

In this paper we primarily focus on the direct proof of 
the heuristic ideas put forward in Refs. [l^Jll|- To do this 
we study in detail the low-energy levels statistics for a 
weakly localized ID Frenkel exciton. The bulk of the pa- 
per is organized as follows. In the next Section the model 
Hamiltonian is described. In Sec. [II the numerical re- 
sults of the low-energy level statistics simulation are dis- 
cussed. Section |^ summarizes the paper and comments 
on the importance of the local discrete energy structure 
for optical response from the disordered Frenkel exciton 
systems. 



II. MODEL HAMILTONIAN 

We consider TV (TV ^ 1) optically active two-level 
molecules forming a regular ID lattice with unity spacing. 
The corresponding Frenkel exciton Hamiltonian reads :E3 
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Here \i) is the state vector of the excited molecule with 
the energy Ei (1 < i,j < Af). The energy is assumed to 
be the Gaussian uncorrected (for different sites) stochas- 
tic variable with zero mean and the standard deviation 
A. On the contrary, the hopping integrals {Ju = 0) 
are considered to be non-fluctuating. These integrals are 
of the dipole-dipole origin: = J/\i — where J is 
the nearest- neighbor coupling. The quantity A/| J| will 
be referred to as the degree of disorder. Hereafter, we 
assume J to be negatiyp, which corresponds to the case 
of the J-aggregates.EB'El In this case the states coupled 
to the light are those at the bottom of the exciton band 
(see, for instance, Ref. ||). 

As it was pointed out in Ref. |l^ the long-range dipole- 
dipole terms strongly affect the exciton-radiative-rate en- 
hancement factor in the presence of disorder. In this pa- 
per we discuss both the nearest-neighbor (NN) approx- 
imation and the effects of the long-range dipole-dipole 
(DD) coupling. 

Let us briefly remind the reader the most important 
features of the one-exciton energy spectrum and eigen- 
functions in the absence of disorder (A = 0). They can 
be found as the solutions of the following eigenvalue prob- 
lem: 

AT 

\K) =Y,V^K^\^) , (2) 

1=1 
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Y,H'WKJ=EK^K^, (3) 

i=l 

where H''^ = {i\H\j), ipK and Ek are the eigenfunction 
and eigenenergy of the one-exciton state \K) respectively; 
the quantity K — ■Kk/{N + 1), where k £ [l,^], plays 
the role of the exciton wavenumber. 

The eigenenergies Ek for the |-e}cact dipole-dipole 
model for A/" S> 1 can be written as:llil 

i^. = -2|J|E^+0(AA-). (4) 

i=l 

The energy spectrum in the NN model can be obtained 
from Eq. H by keeping only the i = 1 term and appears 
to be very different from the exact DD spectrum (see 
below). In contrast to this, theu-eigenfunctions ipK are 
almost the same in both modclsO 

/ 2 

As it was mentioned above we are especially interested 
in the exciton spectrum at the bottom of the exciton 
band, i. e. for if <C 1. In this limit the spectrum was 
obtained and discussed in Ref. |l^: 

Ek =-2\J\Ci3) + \J\l^~\nK]K\ if « 1 , (6) 



where C(3) = X^i^i ~ 1.202. The corresponding ex- 
pression in the NN approximation (the term with i = 1 
in the sum in Eq. (|)) reads: 

Ek = -2\J\ - \J\K^ . (7) 

Straightforward comparison of the two equations shows 
that the long-range DD terms affect the exciton low- 
energy spectrum considerably. First, the exciton band 
bottom is red-shifted in the DD model with respect to the 
NN case: from — 2|J| to approximately — 2.404|J|. Sec- 
ond, the if-dependence of the energy is stronger in the 
DD case due to the logarithmic factor in (^. In particu- 
lar, the energy difference between the two lowest exciton 
states (with K2 = 2n/{Af + 1) and Ki = Tr/{Af + 1)) is 
greater when the long-range terms are considered. 

As well as the eigenfunctions (pK, the oscillator 
strengths are almost the same in both models. For the 
transitions from the ground to lower exciton states they 
can be written as (see, for instance, Ref. ||) 

/a-(Ev-) - '-^i^^ ^=1'2,... (8) 

Here the dipole momentum of the optical transition is 
taken to be unity. Thus, the oscillator strength of the 
transition from the ground state to the lowest exciton 
state {K = Ki = n/{Af+l)), 

/ki = ^(AA+1)«0.81(AA+1), (9) 

is proportional to the number of sites in the chain J\f and 
carries 81% of the total oscillator strength (which is equal 
to Af ). 

III. NUMERICAL RESULTS AND DISCUSSION 

In this section statistical properties of the one- 
exciton low-energy eigenfunctions and eigenenergies of 
the Hamiltonian (|^) are analyzed numerically. We focus 
on the distributions of the localization segment length 
and energy spacings between lower states. The disorder 
degree scaling of the first momenta of these distributions 
will also be the subject of our interest. This requires nu- 
merical solution of the eigenvalue problem (||) for a large 
number of disorder realizations. 



A. Heuristic arguments 

Before doing the detailed statistical analyses let us 
look at a typical realization of the eigenstates (for a 
fixed stochastic distribution of site energies). Figure |l| 
presents an example showing the nine lowest eigenfunc- 
tions ifi, (these are real in our case) and the correspond- 
ing eigenenergies E^,. This set of wave functions was 
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FIG. 1: The typical realization of the low-energy eigenfunc- 
tions of the Hamiltonian (|l|) for the 250 sites linear and the 
disorder degree A/|J| = 0.1. Each plot also shows the cor- 
responding energy of the state Ei, (in units of |J|) and its 
oscillator strength 



obtained for 250 sites chain in the NN model for the dis- 
order degree A/|J| = 0.1. Each plot also displays the 

corresponding oscillator strength f^, ~ ^uij ■ 

Several typical features of the lower eigenstates can 
clearly be seen from Fig. |l|. First, all eigenfunctions are 
well localized at different segments of the chain. The 



lengths of these segments N are smaller than the total 
chain length Af. Some of the eigenfunctions are local- 
ized at the same segment (e. g. manifold of states 4, 6 
and 9) and weakly overlap with the other states. Sec- 
ond, one of the states of each particular (local) manifold 
has no nodes within its localization segment and hence 
have high oscillator strength (e. g. state 4) while the 
other wave functions from the manifold (e. g. states 6 
and 9) have nodes within the segment and low oscilla- 
tor strengths. This allows to classify such states as local 
ground and local excited states respectively. These states 
are similar to the states of a regular non-disordered chain 
(see Eq. (^) at fc = 1 and k = 2, 3). Finally, it is impor- 
tant to note that a local ground state can have higher 
energy than a local excited state which is localized at a 
different segment. For example, the local ground state 
7 has higher energy {Ej — — 2.015| J|) than the local ex- 
cited states 6 {Eq = -2.021| J|) and 5 (£^5 = -2.017| J|). 
This raises the question about the correlations between 
the local ground and excited states from one manifold, 
on the one hand, and between the states from different 
manifolds, on the other hand. 

It is clear from these arguments that grouping the 
states together into local manifolds is a principal prob- 
lem; correct criteria are needed to select the states which 
are localized at the chain segments. When such group- 
ing is done the statistical analysis we are interested in 
becomes simple. 



B. Selection criteria 

To select the states which are localized at the same 
segment first we need to find the local ground states. 
Such states ipi, satisfy the inequality: 



>Co 



(10) 



For the constant Co we use the value Co = 0.95, seeking, 
in other words, for the states which contain at least 95% 
of the wave function density in the main peak. The states 
1 — 4 and 7 in the Fig. |l| meet this criterion and can be 
considered as the local ground states. 

Then we find all local excited states. In order to do 
this we proceed as follows: for each local ground state v 
{ly = 1 — 4, 7 in our example) we seek for the states v' ^ v 
which are localized at the same segment as the state v. 
These (local excited) states should simultaneously meet 
two criteria. The first is similar to the criterion given by 
Eq. (0): 



\vv'' 



>Ci , 



(11) 



where i runs over all sites of the chain. The value for the 
cut-off constant C\ will be discussed later. The second 
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criterion is to be applied to the Fourier transforms 
corresponding to the wave functions ip^,: 



El 

9 



> Co 



(12) 



The later is necessary to drop rapidly oscillating states 
u' , envelope functions of which are localized within the 
same segment as the considered ground state u. The en- 
ergy of such states is in the vicinity of the exciton band 
maximum. Their spatial extension is of the same order 
of magnitude as the localization length of the low-energy 
stateSjSO these high-energy states often meet the crite- 
rion (y). 

The reasonable values for the cut-off parameters Ci 
and 6*2 can be estimated by considering the states of 
a rectangular quantum well. In the limiting case of an 
infinitely deep quantum well the overlap integral for the 
ground state (given by Eq. (||) at fc = 1) and the module 
of the first excited state (Eq. (^) at fc = 2) is equal to 
Im = 8/37r « 0.85. In the more general case of a finite 
rectangular well the value of this integral depends on 
the well depth and width. It varies from about Im (for 
a very deep excited state) to zero. Small values of the 
integral occur for very extended excited states that are 
just appeared in the well. 

Consider the well for which the probability to find the 
particle in the well equals 0.95 for the ground state, then 
for the first excited state this probability equals about 
0.75. The overlap integral of the type (|ll) for these 
states is about 0.8. This value gives an estimate for the 
parameter Ci . The overlap integral of the moduli of the 
Fourier transforms of these states is about 0.56, which 
gives an estimate for the C2. For our estimations we 
considered the limiting case of a rectangular well poten- 
tial which varies infinitely rapidly at the well boundaries. 
In a smoother potential the values of the overlap integrals 
can be smaller. To account for this we use more relaxed 
cut-off criteria with Ci = 0.7 and C2 = 0.5. It is also 
worth noting that straightforward calculations revealed 
only weak quantitative dependence of the final results on 
the values of the cut-off parameters (see Appendix). 



Applying the criteria (|ll|) and (|12|) to our sample set 
of eigenfunctions we find the following pairs of the local 
ground and first excited states: {3,5} and {4,6}. This 
finding looks reasonable (see Fig |l|). The states 4 and 
9 also meet the above criteria. Fig. |l] shows that the 
state 9 has two well defined nodes within the segment of 
localization of the local ground state 4. It therefore can 
be referred to as the second local excited state. However, 
the second and higher local excited states will not be the 
subject of the present study. Such higher excited states 
are usually localized at more extended segments and can 
not be associated with a particular local ground state. 



C. Statistics of the local ground states 

Figure |l| clearly shows that the spatial extension of 
a local ground state (its localization length N) fluc- 
tuates from one state to another. Different quan- 
tities can be used to characterize the extension of 
a wave function: (i) the inverse participation ratio 
(IPR), defined as |<(5i/i|'', (ii) the mean-square dis- 

placement |v3^^p(i - i)^] , where i = Y.i 
(iii) the inverse of the Lyapunov -exponent X,y — 
— limi_»oo (l/2i) In |(y9i,ip and others.cl Although these 
quantities yield slightly different localization lengths, 
they are of the order of the number of sites over which 
the wave function is extended. 

For the local ground state, an alternative quantity can 
be used as the measure of its extension. According to 
Eq. (||) , the oscillator strength of the ground state of the 
ideal linear chain is about the same as the chain length. 
In other words, /i/ also carries information about the ex- 
tension of the corresponding local ground state wavefunc- 
tion (p^. Furthermore, the oscillator strength can be mea- 
sured experimentally, e.g. it can bej-determined from the 
exciton spontaneous decay kineticsoEj) . For this reason 
we use the oscillator strength as the measure of localiza- 
tion of the local ground states. 

The number of sites in the chain M should be much 
larger than the typical oscillator strength /* (M ^ /*). 
As /* was found to be of the order of several tens for 
the disorder degree A/|J| = 0.1 and decreases with in- 
creasing A/I J|, we consider 400 cites chains (M — 400). 
Applying the selection criterion ( |l0| ) we calculated the 
/-distribution for the local ground states. The result ob- 
tained in the NN approximation is presented in the Fig. |^. 
Figure |^( a) shows the (typical) oscillator strength distri- 
bution P(/) obtained for a 400 sites chain by averaging 
over 5000 realizations of disorder for the disorder degree 
A/I J| — 0.1. The main features of the distribution are 
clearly seen in the figure: (i) - asymmetry with respect 
to the expectation value of /, (ii) - long tail in the large 
lengths region, and, on the contrary, (iii) - steep drop off 
in the region of small values of /. 

Figure ||(6) presents the disorder degree scaling of 
the first momentum of the -P(/) distribution, /* = 
/ dfP{f)f. Error bars in Fig. ||(6) show the standard de- 
viation CT EE [/ dfP{f){f - /*)2]i/2. The best fit (dashed 
line) to the numerical data (diamonds) is given by the 
power-law approximation in the form 



-0.67 



(13) 



The theoretical estimate for the localization length was 
presented in Refs. |l^,|ll| (see also Ref. ^ and was based 
on the following arguments: An exciton localized at a 
segment of size N* feels_a reduced (exchange-narrowed) 
disorder a = A/\/iV*jia this reduced disorder degree 
must be of the same order of magnitude as the energy 
separation between the two lowest local exciton levels. 
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FIG. 2: (a) - The distribution P{f) of the local ground state 
oscillator strength / obtained for a 400 sites chain and the dis- 
order degree A/|J| = 0.1 (5000 realizations of disorder was 
considered, each point of the distribution is averaged over the 
interval (accumulation length) 5 = 1). {b) - the disorder de- 
gree scaling of the first momentum f* = J dfP{f)f obtained 
for a 400 sites chain (300 realizations of disorder was consid- 
ered for each value of the disorder degree A/| J|). Diamonds 
present the numerical data while the dashed line gives the 
best power-law fit. Error bars show the standard deviation of 
the distribution, a = [J dfP{f){f - ffY''^. 



El is smaller than 



El -El = 37r2|J|/7V*. Indeed, if £;* 
a the two levels would be strongly mixed by disorder and 
the corresponding exciton wave functions would tend to 
reduce their extension. On the contrary, if E^ — El ex- 
ceeds a then the disorder is perturbative and the wave 
functions would tend to increase their extension. Thus, 
the equality E^ — El = a gives the self-consistent re- 
lationship between the spatial extension of the exciton 
wave functions and the disorder degree which the local- 
ized exciton feels. This relationship provides an estimate 
for the localization length N*: 



N* 



3^VI 
A 



2/3 



9.57 
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Ml 



-2/3 



(14) 



This theoretical scaling law is in good agreement with the 
scaling law of /* obtained numerically, the latter being 
smaller by the factor of about 2. This is not surprising 
because, according to Eq. (||), the oscillator strength is 
smaller than the localization length. The mean localiza- 
tion length N* was also calculated {N was defined as 



the number of sites under the main peak of the wave- 
function, for which |(^yi|/maxi{|(/3!yi|} > 0.1). The mean 
oscillator strength /* was found to be proportional to 
N*: f* = cN*, where c « 0.66 . Rewriting Eq. (|l|) in 
terms of N* yields TV* « 8.71 (A/| J|)-0-67, which is in 
excellent agreement with the theoretical estimate (p^. 



D. Statistics of the energy levels near the lower 
exciton band edge 

In this section we analyze two eigenenergy distribu- 
tions which are of primary interest for us. The first is 
the distribution of the energy spacing between the local 
ground states of different segments, while the second is 
the distribution of the energy spacing between the lo- 
cal excited and ground states of the same segment. We 
denote them as Pii{SE) and P2ii SE) re spectivel y . Uti- 
lizing the selection criteria of Sec. |IIIB (Eqs. (|To|)-(|T2|)) 
we calculated these distributions. The results of these 
calculations are presented in Figs. |3| and ^. 

Figures ||(a) and ||(&) present the (typical) distribu- 
tion of the energy spacing between the local ground 
states Pii{5E) obtained for a 400 sites chain by aver- 
aging over 5000 realization of disorder. The disorder de- 
gree A/I J| — 0.1 and the NN model was used. The plot 
(a) shows Pii{SE) for the adjacent (most overlapping) 
ground states only. These states were selected by the cri- 
terion l^vil > 0.01. The distribution presented 
in Fig. |3|(&), was obtained without the above restriction, 
i. e. for all ground states. 

The two distributions are very close to each other. The 
most important feature of both is that they are nearly 
symmetric with respect to zero energy spacing, in other 
words, they are characterized by zero mean. The disorder 
degree scaling of the corresponding first momenta (not 
presented here) confirms this fact. A weakly pronounced 
downfall in the vicinity of the zero energy spacing of the 
distribution for the adjacent states (plot (a)) is due to the 
repulsion of close energy levels of the overlapping states. 
Similarly, two degenerate states in a camel's back po- 
tential are split into a doublet when interaction between 
them is considered. This downfall tends to disappear 
when all ground state are considered (see Fig. ||(6)). 

The plot (c) of Fig. ^ shows the disorder degree 
scaling of the standard deviation of the distribution 
Piii6E), ail = [J d{SE){SEfPii{SE)]^/'^, calculated for 
all ground states. The best power-law fit (dashed line) 
to the numerical data (diamonds) is given by 



0-11 « 0.67 I J| 
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Ml 



1.33 



(15) 



The scaling law of an is close to-iher-scaling law of the 
one-exciton absorption linewidthE]'E2l'E3, which supports 
the idea that these quantities are proportional to each 
other. 
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FIG. 3: {a,b) - The distributions of the energy spacing be- 
tween the local ground states obtained for a 400 sites chain 
and the disorder degree A/| J| = 0.1 (5000 reahzations of dis- 
order was considered, each point of the distributions is aver- 
aged over the interval (accumulation length) 5 = 0.0015). The 
plot (a) shows the distribution Pii(5E) for the adjacent lo- 
cal ground states (with the overlap integral \ > 
0.01). The distribution shown in the plot (fe) was calculated 
for all local ground states, (c) - the disorder degree scaling of 
the standard deviation an = [f d{SE){SEfPii{SE)]^^^ ob- 
tained for a 400 sites chain (300 realizations of disorder was 
considered for each point). Diamonds present the numerical 
data while the dashed line gives the best power-law fit. All 
energies are in units of |Jj. 



Figure |4|( a) shows the (typical) distribution P2i{6E) 
of the energy spacing between the ground and first ex- 
cited states localized within the same segment. This dis- 
tribution was obtained for the same set of parameters 
[M = 400, A/I J| = 0.1, 10000 realizations). In contrast 
to the Pii{SE) (Fig. |(a, b)) the distribution P2i{SE) is 



asymmetric with respect to the expectation value of SE. 
The high-energy-spacing shoulder of this distribution is 
much longer than the shoulder in the region of small val- 
ues of SE. However, the most important features of the 
P2i{SE) are the following. First, it is characterized by 
non-zero mean, which unambiguously shows that the en- 
ergies of the pairs of states that are localized within the 
same segment are highly correlated, in the sense that 
they undergo repulsion resulting in the local discrete en- 
ergy structure. Second, the distribution P2i{5E) drops 
down rapidly at 5E ~ 5Ec and vanishes below the en- 
ergy 5E(.. For the particular value of the disorder degree 
A/I J| = 0.1 this energy is small compared to the char- 
acteristic energy scale in the figure {5Ec ~ 0.001). The 
log-linear insert presents a blow-up of the distribution 
function P2i{5E) in the region of very small energy spac- 
ings 5E ^ SEc (each point is averaged over the interval 
S = 0.0005 in the insert). 

The quantity of interest in the present case is the 
first momentum of the distribution P2i{5E). Denote it 
as E2 — E^. The physical meaning of this quantity is 
clear: it is the mean energy spacing between the levels 
of a local manifold. The calculated disorder scaling of 
E* ~ E* ^ J d{6E)P2i{SE)6E (diamonds) is shown in 
Fig. ^(&) together with the best power-law fit (dashed 
line) to the numerical data: 



E*2 - El 



0.40 |J| ( ^ 



1.36 



(16) 



Note that the disorder scaling laws for E2 — El and cth 
are very close to each other. The theoretical estimate of 



4/3 



the E2 — El scaling was obtained in Refs. |10|,|11| 



e;-ei 



4/3 



0.32 |J| ( A 



(17) 

and agrees excellently with the calculated behavior. This 
result was also obtained numerically in Ref. |lj, where 
the detuning between bleaching and induced absorption 
of the pump-probe spectrum was simulated. This de- 
tunin g i a a good measure for the exciton delocalization 
lengthEiy. 

It is remarkable that rescaling of the numerical data 
for El - El in terms of /* gives E^ - El - (l//*)^ 
the proportionality that holds for an ideal linear chain 
of length /*. This is yet another confirmation of the 
existence of the local discrete energy structure of a weakly 
localized Frenkel exciton near the lower band edge. 

It is also of interest to compare the mean E2 — El of 
the P2i{5E) distribution with the standard deviation an 
of the Pii{SE) distribution. The figures ||( c) and^(&) as 
well as the equations ( p^ and (17) for these quantities 
show that E2 — El < ail. This finding clarifies the fact 
that the local discrete energy structure can not be seen 
neither in the density of states (see Refs. [7| jr9|,^,p3|) nor 
in the linear absorption spectra (see Refs. [7|, |13yi5| , |19| ); the 
energy fluctuations of the local ground states exceed the 
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FIG. 4: (a) - The distribution of the energy spacing be- 
tween the local ground and first excited states of the same 
segment, P2i{5E), obtained for a 400 sites chain and the dis- 
order degree A/|J| — 0.1 (10000 reaUzations of disorder is 
considered, each point of the distributions is averaged over 
the interval (accumulation length) 5 — 0.0025). The insert 
shows the behavior of the distribution function P2i{SE) in 
the region of very small energy spacings SE (each point is av- 
eraged over the interval S = 0.0005) (fe) - The disorder degree 
scaling of the first momentum E2 — El = J d{5E)P2i{5E)5E 
obtained for a 400 sites chain (each point is averaged over 
300 reahzations of disorder). Diamonds present the numer- 
ical data while the dashed line gives the best power-law fit. 
Error bars show the standard deviation of the distribution, 
(721 = [/ d{5E)P2i{5E){E2 - El - SEf]^/"^. AH energies are 
in units of \J\. 



characteristic energy scale of the local energy structure, 
hence, the fine local structure is hidden by the inhomo- 
geneous broadening of the local ground states. 

The effect of level repulsion is well known from the. 
studies of spectra of complex systems, such as nuclei. O 
These studies, are often based on the random matrix the- 
ory (RMT)o In particular, this theory predicts level re- 
pulsion for the eigenvalues of a real symmetric matrix 
chosen at random from the Gaussian orthogonal ensem- 
ble, i.e. when all matrix elements are Gaussian stochastic 
variables statistically independent of each other. Despite 
a seeming similarity of this result and ours they cannot 



be directly related to each other. First, the off-diagonal 
matrix element of the Hamiltonian (|l|) are not stochas- 
tic variables while they are in the random matrix theory. 
Moreover, in the RMT the off-diagonal elements have 
zero mean, while in the Hamiltonian (^ they are all neg- 
ative. Second, we found the well-pronounced effect of 
level repulsion only for a sub-ensemble of levels. Simu- 
lations for the entire set of levels show almost no level 
repulsion in the vicinity of the band bottom (see Fig. |^). 



E. Long-range coupling effects 

Complete set of the above mentioned simulations was 
also performed for the exact DD coupling model. It was 
found that /*, cth, and E2 — El are strongly affected 
by the long-range dipole-dipole terms. The following ex- 
pressions give the best fit to the numerical data: 



/* « 9.74 



cTii « 0.50 |J| 



A 



1.36 



E* -El ^0.2,1 \J\ ( A 



1.53 



(18) 



(19) 



(20) 



Unlike the NN model case, the rescaling formula E2 — 
El r-. (l//*)2 does not hold for the exact DD model. This 
finding is not surprising, as it agrees with the change in 
the exciton energy spectrum of a regular chain in the DD 
approximation as compared to that obtained in the NN 
approximation (see Eqs. (^ and (0)). The logarithmic 
factor in Eq. (|) is breaking the E* - El ~ (1//*)^ de- 
pendence. The deviation of the exponents 0.80 and 1.53 
in Eqs. (|l|) and (|^) from theJ^N model values 2/3 and 
4/3 is also due to this factorJld As a matter of fact, the 
power-law approximation is not adequate for the exact 
model. Using the exact low-energy spectrum of the ID 
localized exciton and the rule described at the end of 
Sec. IIIC), the authors of Ref. |l^ obtained the correct 
equation for the localization length N*: 



iV*3/2 

In iV* + ( 



(21) 



with a = 3/2 - ln{2^/^TT) w -0.57 and h = in'^ 30. 

Straightforward comparison of the scaling laws of the 
quantities /* and N* (the latter is defined as the number 
of sites under t he ma in peak of the wavefunction, see also 
the end of Sec. [11 Cj ), shows that, like in the NN model, 
these quantities are proportional to each other f*—c N* , 
where c « 0.62. Fitting Eq. to the numerical data 
for f* , the following optimal values of the parameters 
were found: a = —0.54 and b — 25.1. These values are in 
excellent agreement with the theory. 
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FIG. 5: The distribution of the energy spacing between levels 
in the vicinity of the exciton band bottom, Pij{5E), obtained 
for a 400 sites chain and the disorder degree A/|J| = 0.1 
(5000 realizations of disorder was considered, each point of 
the distributions is averaged over the interval (accumulation 
length) 5 — 0.00015). All energies are in units of |J|. Note 
that that this distribution shows almost no level repulsion. 



In conclusion to this section we would like to com- 
ment on the disorder scaling law of the exciton-radiative- 
rate enhancement factor which was obtained numerically 
in Ref. ^ for the exact dipole-dipole coupling model. 
This factor is introduced as follows: define the aver- 
age oscillator strength per state at energy E: n'^{E) = 
A{E)/p{E), where A{E) and p{E) are the absorption 
spectrum and the density of exciton states respectively, 
then maxlfj,'^ (E)] is an effective measure for the enhance- 
ment of the radiative rate. This quantity also contains 
information about the spjaiial extension of the exciton 
states coupled to the light .Ej The authors of Ref. |l^ found 
from their simulations that max[fi^{E)] ~ (| J|/A)"-^''. 
The exponent of this scaling law, 0.74, differs from the 
value obtained in our simulations, 0.80 (Eq. (ITsI)). We 
believe that this discrepancy of the order of 10% origi- 
nates in the different definitions of /* used in our study 
and in Ref. ITgl. 



N* (localization length). The lowest state of each set has 
no nodes within the localization segment and therefore 
can be interpreted as the local ground state, while the 
next state in a set has a node within the segment and 
is the first local excited state, etc. The average energy 
spacing of the local ground and first excited states ob- 
tained within the framework of the NN model follows 
the (1/iV*)^ dependence that holds for an ideal chain of 
length N*. 

On the contrary, the energy spacing of the non- 
overlapping local ground states is characterized by zero 
mean and distributed in the energy interval that is of 
the same order of magnitude as the mean energy spac- 
ing of the local energy structure. Therefore, the lo- 
cal energy stj.iir.ture appears to be hidden in the den- 
sity Qf-stetefiLrlljCj and in linear absorption spectra as 
well.l3ll3ll^£3 Nevertheless, this structure can be revealed 
by the-|pump-probe spectroscopy sharing two-exciton 
states.Oli Jijjdfied, due to fermionic nature of ID Frenkel 
excitons,L§E3cZl'L3 one-to-two and zero-to-one exciton op- 
tical transitions involving the states of the same local 
structure would be blue-shifted with respect to each 
other. This shift is equal to the energy spacing between 
the local ground and first excited states. The blue shift of 
one-to-two and zero-to-one exciton optical transitions in 
J-aggregates ot^jseudoisocyanine bromide was observed 
experimentallylljO providing unambiguous confirmation 
of the existence of the local discrete energy structure in 
the vicinity of the exciton band minimum. 
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IV. SUMMARY AND CONCLUDING 
REMARKS 



APPENDIX A: EFFECT OF THE CUT-OFF 
PARAMETERS 



In this paper we discuss the statistics of energy levels 
of the ID Frenkel Hamiltonian with on-site randomness 
near the lower band edge. The heuristic arguments on 
the existence of the local discrete energy structure, which 
were put forward in Refs. |l^ and ^are confirmed by the 
detailed statistical study. Selecting the states by means 
of the overlap integrals we find that the energy levels of 
the well overlapping states indeed undergo repulsion re- 
sulting in the local energy structure which is similar to 
the structure of an ideal linear chain of the reduced length 



In this section the weak dependence of the results on 
the cut-off parameters C\ and C2 is illustrated. Figure ^ 
shows several distributions P2i{SE) calculated for differ- 
ent sets of the cut-off parameters: Ci — 0.7, C2 = 0.3 
— dash-dotted line, Ci = 0.7, C2 — 0.5 — solid line, 
Ci = 0.8, C2 = 0.3 — dotted line, Ci = 0.8, C2 = 0.5 — 
dashed line. These distributions are very close to each 
other, which results in only weak (less than 15%) quanti- 
tative dependence of distribution momenta on the values 
of the cut-off parameters. 
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FIG. 6: The distribution of the energy spacing between the 
local ground and first excited states of the same segment, 
P21 (SE), obtained for a 200 sites chain and the disorder degree 
A/|J| = 0.1 (40000 realizations of disorder was considered, 
each point of the distributions is averaged over the interval 
6 = 0.0025). The curves are calculated for different sets of 
the cut-off parameters: Ci = 0.7, C2 = 0.3 — dash-dotted 
line, Ci = 0.7, C2 = 0.5 — solid line, Ci = 0.8, C2 = 0.3 — 
dotted line, Ci — 0.8, C2 = 0.5 — dashed line. All energies 
are in units of |J|. 
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